Spin swap vs. double occupancy in quantum gates 
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We propose an approach to realize quantum gates with electron spins localized in a semiconductor 
that uses double occupancy to advantage. With a fast (non-adiabatic) time control of the tunnelling, 
the probability of double occupancy is first increased and then brought back exactly to zero. The 
quantum phase built in this process can be exploited to realize fast quantum operations. We 
illustrate the idea focusing on the half-swap operation, which is the key two-qubit operation needed 
to build a CNOT gate. 

PACS numbers: 75.10.Pq 
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Many solid state proposals for quantum computation 
make use of the spin of an electron localized in a quan- 
tum dot 0,13] or by a donor ion 3| as a qubit. In these 
schemes the overlap between the electron wave functions 
in neighboring dots is controlled by gates, and the ex- 
change interaction provides the mechanism for the real- 
ization of conditional operations on two qubits. Among 
the requirements for such proposals is that one wants no 
or small double occupancy, i.e. small probability to have 
the two electrons in the same dot or donor. In fact, the 
state space S for 2 qubits has to be 4-dimensional and 
the mixing with the states of double occupancy brings 
the physical system out of the computational space. Hu 
and Das Sarma jj] observed that the double occupancy 
will in principle degrade any implementation of quan- 
tum gates, and noted the intrinsic conflict of trying to 
minimize this double occupancy while making the over- 
lap large enough to realize a sizeable two-qubit opera- 
tion, like e.g. a spin swap. The problem was further 
analyzed in Ref. [a], where it was shown that using an 
adiabatic control in the quantum dots case, the double 
occupancy was small enough so that the degradation can 
be made negligible, consistently with the requirement of 
large enough overlap. The operation critical to quan- 
tum computing is the half-swap operation (also written 
as V SWAP), i.e. having the state going half of the way 
to the complete swap, and we focus on this process to 
illustrate the approach. 

We show here that a half-swap operation without degra- 
dation can be realized by designing a time control where 
the probability of double occupancy first increases and 
then is brought back exactly to zero. In this cycle the 
quantum state in S picks up a phase that is exploited 
in completing the process. As a consequence, the swap 
can be much faster than that in the case where double 
occupancy is kept small at all times. We will consider ex- 
plicitly two cases where the control of the tunneling am- 
plitude follows in time a square and a hyperbolic secant 
pulse shape. The conditions for a perfect v SWAP, i. e. 
a fast half-swap gate with double occupancy identically 
zero, can be given in these two cases by simple analytic 



rules. The results are exact within the single band model 
considered, the validity of which will be discussed below 
for realistic systems. 

The dynamics during the half swap can be described 
by using two spatial orbitals, which we call w a and Wb 
(w for Wannier functions), orthonormal, and localized 
respectively around the two sites a and b 0, 0, |(| . The 
w's are linear combinations of non-orthogonal "atomic" 
functions whose overlap is A. This implies a 4-state one- 
electron basis, and therefore a 6-dimensional complete 
vector space for 2 electrons, pj The Hamiltonian H con- 
tains the kinetic and potential energies, plus a magnetic 
field interacting with the orbital magnetic moment. The 
quantum control is realized by changing the profile of the 
potential energy, therefore modifying the tunnelling am- 
plitude between the two sites. We will not consider a 
magnetic field term, since it is irrelevant in the proposed 
scheme. The 2-electron basis used is 
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where c\ Ja creates an electron in the orbital w v with 
spin a. Physically we speak of <&„ with v = 1,4,5,6 
as having 1 electron near each site, while for v = 2,3, 
the states have 2 electrons on one site, and are called 
doubly-occupied states. Assuming inversion symmetry, 
the time-dependent Hamiltonian Hij(t) = ($>i\H(t)\$>j) 
reduces by symmetry to four 1 x l's (the triplet i/44, H55, 
and H 66 , and the ungerade singlet H 33 ), plus a 2 x 2 
which yields the remaining two (gerade) singlets. Wc 
assume that before applying the control pulse, at time 
t = —00, the wave function is |\&(— 00)) = c£*c£i|0) = 
2~ 1/2 (|$i) + |*5)) • Given the fo rm of the initial state 
and the symmetry properties discussed above, ^3), ($4) 
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and | $6) do not enter in the dynamics and we calcu- 
late =a(t)\$i) + b(t)\$ 2 )+c(t)\$5) by the (time- 
dependent) Schrodinger equation 
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where M2(i) = 2H 12 (t), hA = H 22 - H n and hJ P (t) = 
Htt(t) — Hn(t). Notice that tt(t) is first order in the 
overlap A, while A depends only on the on-site and long- 
range Coulomb interaction: hA = (w a , w a \v\w a , w a ) — 
(w a , wi,\v\w a , Wb), where v is the Coulomb interaction 
between the two electrons. This term is of order A , 
and therefore only weakly dependent on the tunnelling 
amplitude. So we assume A constant during the gate 
operation. hJpjt) — — 2(w a , Wb\v\wb, w a ) is the poten- 
tial exchange y|, which is 0(A 2 ). The hopping term 
htt(t) = A{w a \h\ \ wb)+4:(wa, Wb\v\w a , w a ), where hi is the 
single particle kinetic and potential energy, is responsible 
for the kinetic exchange. Both Jp(t) and £l(t) are con- 
trolled by the tunnelling amplitude through their depen- 
dence on the overlap A. We remark that these parameters 
simplify in the Hubbard model to Jp(t) — and hA = U. 
In this case the hopping amplitude £l(t) is the only mech- 
anism for the exchange interaction (Jk — — h 2 Vl{t) 2 /2U 
for small A, the kinetic exchange [8J). The presence of 
the potential exchange Jp{t) is not essential in the real- 
ization of the proposed scheme, but we keep it in order 
to prove the reliability of the scheme in the most general 
case. 

It is now simple to calculate the quantities we want to 
monitor. One is the expectation value of s z — n a -\ — n a i, 
the z-component of spin at site a. (n c are occupation 
numbers.) It is sufficient to consider just one site since 
the total z-component of spin is conserved and so must 
remain 0. Using the property s z |$i) = |^s), s z \&5) = 
|$i), and s z |$ 2 ) = 0, we find 



a z (t) = (*(t)|s,|*(t)) = 2Re[a(t)c*(t)} 



(3) 



We also want the probability Pd of double occupancy, 
i.e. the probability of finding the system in the state $2. 
This is 



P d {t) = |<$ 2 |*(t))| 2 = |6(t)| s 



(4) 



In order to have the perfect \J SWAP at the end of the 
pulse, i.e. at t = +00, we need to design the control in 
such a way that the two conditions 



P d (+cx>) = 
o-4+oo) = 



(5a) 
(5b) 



are verified. One can see that Eq. © implies |*(+oo)) 



2~ 1/2 (<4 r 4 1 ± i4j_4 T )|Q) (Bell entangled states). In 
the usual approach a weaker constraint than the one in 



Eq. (|5a|l is used (i.e. Pd(+oo) <§; 1) satisfied by slow 
adiabatic switching and small mixing with doubly occu- 
pied states. This condition is implicit when we describe 
the system of two localized electrons using the (effec- 
tive) Heisenberg Hamiltonian. Our approach consists in 
finding first the general conditions on the control to sat- 
isfy exactly Eq. (|5al) . which does not imply necessarily 
adiabaticity and small average mixing with the doubly 
occupied states, and then adjust the remaining control 
parameters to satisfy Eq. I|5bjl . 

The time dependence of c(t) can be integrated directly 
from the Schroedinger equation and we obtain 



c(+oo) 
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The dynamics of a(t) and b(t) can be mapped to the evo- 
lution of a pseudospin in the presence of a time dependent 
effective magnetic field. The condition in Eq. I|5al) means 
that the pseudospin makes a closed loop in the Bloch 
sphere, while the condition in Eq. (|5b(l ensures that the 
phase picked up by the pseudospin evolution during such 
a closed loop compensates the phase due to the potential 
exchange for a perfect V SWAP. The analytic solution 
for the evolution of a spin in a time dependent magnetic 
field can be found only in particular cases , and we will 
discuss here two of them: the Rabi solution 0|j corre- 
sponding to a square pulse shape, and the Rosen and 
Zener solution corresponding to a hyperbolic secant 
pulse shape. 

The Rabi solution describes the precession of the pseu- 
dospin in a constant effective magnetic field. A square 
pulse in the tunnelling amplitude keeps this effective field 
on for a time T, and then turns it off. We can write the 
hopping term as a constant Q(t) — O, and consider the 
time evolution given by the Rabi solution from a time 
t = to t = T (just the solution for t-independent H). 
In a similar way the potential exchange term Jp(t) = Jp 
for < t < T and zero otherwise. The probability of 
double occupancy at the end of the pulse is given by: 



P d (T) = \b(T)f = 



n 2 /2 



() 2 
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Vn 2 + 5* 



(7) 



Therefore the condition in Eq. H5a() is satisfied by fixing 
the length of the pulse to 



Tn 



2mr 



Vft 2 + A 2 



(8) 



The integer n identifies the number of full precessions 
of the pseudospin during the pulse. We find from the 
Rabi solution that a(T n ) = exp(-iAT„/2)/\/2. Using 
this solution in Eq. ©, the condition in Eq. J5b|) for the 
perfect V SWAP is rewritten as 



( A - 2 J P 
cos mr = = 0. 

VV^ 2 + A 2 



(9) 
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FIG. 1: cr z (t) and Pd(t) for two electrons confined by hydro- 
genic potentials. Infinite tunnelling barriers are switched off 
during the time interval < t < T. Time t is in units of 
h/Ky* . (a) The separation between the sites is 4.53 au*. At 
the end of the pulse we have a half swap with a small but 
nonzero Pd- (b) The separation is 2.38 au* and the energy 
levels satisfy the condition in Eq. IIUI . The half swap is real- 
ized in a shorter time and with Pd = identically at the end 
of the pulse (see inset). 



This can be recast in terms of the separation between 
the triplet and the lower singlet eigenvalues J, and the 
separation between the two gerade singlets AE as 



2nJ 
~AE 



2k + 1 



(10) 



where k is any integer. 

To illustrate the idea we show in Fig. ^ the evolution 
of a(t) and Pd(t) in the square pulse case calculated as- 
suming a hydrogenic confinement potential and localized 
functions of the Slater form, which models well the case 
of electrons confined in shallow donors. We want to stress 
that the general approach we are suggesting is indepen- 
dent of this particular choice. In fact, for any reasonable 
form of the confinement potential it is possible to modify 
the ratio in Eq. I|10|l by changing the maximum tunnelling 
amplitude or the separation between the dots/impuritics. 
When the separation is large, so the overlap is small (or 
the tunnelling barrier is large), there is a small admixture 
of doubly-occupied states at any time. This is the region 



where the low-lying states are described accurately by the 
Heisenberg Hamiltonian 2 JS a • S^. The spins here are as- 
sociated with the sites; e.g. 5+ = e^c a j 8}. This case 
is considered in Fig. ^ (a), obtained with a separation 
between the localized electrons of 4.53 au*. We indicate 
by au* and Ry* the atomic units scaled by the electron 
effective mass and the static dielectric constant. Time is 
in units of h/~Ry*. For a typical GaAs based system 1 
Ry* - 5 meV, 1 au* - 80 A, and h/Ry* ~ 0.1 ps. We 
see that at the end of the pulse we satisfy the condition 
for the half swap a z = with a small double occupancy 
Pd ~ 0.03 (see inset in Fig. 1 (a)). The new approach is 
illustrated in Fig. ^ (b) corresponding to a separation of 
2.38 au*. At this separation the energy levels satisfy the 
condition in Eq. (|10(l yielding a perfect half swap with 
both <r 2 = and Pd — identically. Fig. ^ (b) shows a 
marked increase in swap speed and a much higher prob- 
ability of double occupancy during the gate compared to 
Fig. E(a). The case in Fig.Q](b) corresponds to n = 1, 
higher n will give a perfect half swap in longer times. 

The use of square pulses gives a direct illustration of 
the new approach. However, for practical purposes, a 
continuous time dependence of the pulse is desirable. Us- 
ing the picture of the pseudospin in an effective magnetic 
field, we need to design in the general case a B e //(i) that 
generates closed loops in the Bloch sphere. Once this 
first condition is realized, the condition on the control 
parameters for a perfect half swap can be found accord- 
ingly. This is possible for the hopping amplitude of the 
form n(t) = n/ cosh (nt/T). Then Sl{t) dt = fiT, 
the same as for the square pulse. The probability of dou- 
ble occupancy at the end of the pulse (at t = +oo) can 
be obtained from the Rosen and Zener result 



Pd = 



sin 2 nr/2 

2 cosh 2 AT/2' 



(11) 



In the adiabatic approach the condition AT 3> 1 in the 
denominator makes this probability small. Our approach 
takes advantage of the numerator in Eq. (|11|> . which im- 
plies that for 



2nn 



(12) 



we will have no double occupancy at the end of the pulse. 
In order to determine the phase of a(+oo), we use the 
transformation 2z = tanh(7rf/T) -I- 1 to map the infinite 
time interval to [0, 1], and the equation for the coefficient 
a(z) reads 



(z 2 - z)a" + ± + 



iAT 
~2~^~ 



z)a' + 



~2v7 



a = . (13) 



This is the hypergeometric equation and in general we 
obtain a solution in terms of a linear combination of hy- 
pergeometric functions. However, it is interesting to no- 
tice that if the condition for no double occupancy is met, 
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FIG. 2: a z (t) and Pd(t) for two electrons confined by hydro- 
genic potentials with barriers controlled using a hyperbolic 
secant pulse shape. The separation is 2.36- ■ ■ au*. Time t is 
in units of R/Ry*. The inset shows the Log of o~(t) and Pd(t) 
at longer times. 

then (VlT/2'k) = n, and we will fall in the degenerate case 
of the hypergeometric equation. As discussed in Ref. [l^ . 
the solution in this case consists of a polynomial of degree 
n. We find that the solution to Eq. l|13|l with the initial 
condition a(Q) = 1 and with |a(l)| 2 = 1 can be written 

explicitly (z) — n\P n "' "(1 — 2z)/(a n + 1)„ , where 
P^' 1 (x) are Jacobi polynomials, a n = — 1/2 + in A/ft, 
and (x) n = r (x + n) /T (x). The condition in Eq. (|5b|) 
can then be written in a compact form using the hyper- 
geometric function (which in this case terminates after n 
terms in the series) as 

Re[F(-n, n, 1/2 + mA/O, \) e iinJ ^/ n ] = , (14) 

which is the parallel of Eq. (|5J) for this continuous pulse 
shape. The fact that the potential exchange Jp(t) is 
0(A 2 ) led us to take J P (t) = J P sech 2 (t/T) in Eq. <|T1). 
although the solvability and essential aspects of the so- 
lution are not limited to this. We show in Fig. the 
realization of a perfect half swap in the hydrogenic case. 
At the separation of about 2.36 au* the matrix elements 
of the Hamiltonian satisfy the condition in Eq. I|14|l and a 
perfect half swap can be realized. The case in Fig. [^cor- 
responds to n = 1 and we have checked that for higher 
n we can always find at least one separation to satisfy 
Eq. (|14fl . The inset shows the exponential converging of 
Pd and a z toward the desired value. We remark that in 
GaAs based systems this translates in a half swap oper- 
ated with an accuracy of 10~ 10 in a time of the order of 
3 ps. 

We now consider likely sources of error and other dif- 
ficulties to be expected in a real implementation of the 
ideas proposed. We pose the question, if and A deviate 
slightly from values needed for a perfect half-swap, how 
much do er(+oo) and Pd(+oo) differ from their ideal val- 
ues of zero? Using the analytic solutions we proved that if 



the deviation of the parameters is 0(e) then the deviation 
ct(+oo) =0(e) andPd(+oo) =0(e 2 ). Higher intra-dot ex- 
citations will complicate our considerations. If they are 
large compared to hA, however, their effect will be small. 
To get an estimate of the requirements, consider these ex- 
citations for the infinite square well, side length L. Then 
SE = 3ir 2 h 2 /2m*L 2 = 67r 2 (au*/L) 2 Ry*, is the first exci- 
tation energy, and TiA w 2e 2 /(eL) = 2(au*/L)Ry*, so 
SE/hA w 37r 2 (au*/L). For two GaAs quantum dots 
of lateral size L = lOOA, this gives SE/hA w 27 and 
TiA ss 2Ry*, satisfying the requirement SE/hA >> 1. 
Finally, we remark that in order to take advantage of the 
proposed scheme the time resolution needed for control- 
ling the gate is related to the intra-dot Coulomb repulsion 
hA. In principle, a subpicosecond switching time for a 
pulse gate is not out of the question By increas- 

ing the dot size L while maintaining a large SE/hA it 
is possible to increase the resolution time A -1 . We esti- 
mate that this condition can be well satisfied for large (~ 
300-400 A) GaAs quantum dots in such a way to limit 
the resolution to the picosecond range, which is within 
state-of-the-art capabilities. 

In conclusion, we have shown that by non-adiabatic 
control it is possible to take advantage of the double oc- 
cupied states to make fast and accurate half swap oper- 
ations. The idea can be applied to many systems cur- 
rently under investigation for experimental realization of 
quantum computing. Even in systems where the mixing 
with double occupancy is small, the simple rules given 
here can be used as an error correction method. Possi- 
ble generalizations to more than two spins fl4l | are under 
consideration. It is also of interest to consider more gen- 
eral pulse shapes, particularly with the object of seeing 
what forms of B e y f (t) will allow the pseudospin to return 
to its starting value (i.e. corresponding to Pd = 0). 
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